REPORT  DOCUMENTATION  PAGE 

MjL  1  J*"*  MJL  ^  a  I  estimated  to  swe^a^e  1  hour  per  response 

•  *  mm  m  rewtewino  the  coiienion  of  informati 


form  Approv#</ 
0MB  No.  0704  018B 


estimated  to  swera^e  1  hour  per  response,  ir^dudin^  the  time  tor  reviewing  instructions,  searching  eMsting  data  sources. 
I  and  reviewing  the  coiieciion  of  information  Send  comrnents  regarding  this  burden  estimate  or  any  other  aspect  of  this 
}  this  burden  to  WMhifsgton  Headduarters  Services.  Directorate  for  information  Operations  and  Heports.  t2)S  ietferson 
0  the  Office  of  Management  and  Budget.  Paperwork  Acdurtion  Prfuect  (O704«O  f 8i>.  Washington^  DC  20503- 


URORT  DATE 


4.  TITU  AND  SUBTITLE 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final  Report  lOOct  92  -  30  Sep  93 


5.  FUNDING  NUMBERS 


Computer  simulations  of  radiation  generation  from 
relativistic  electron  beans 


ntiiiiiwnoi 


Dr  Anthony  T.  Lin 


7.  PERFORMING  ORGANIZATION  NAMECSTANO  A00RESS<ES) 
University  of  California 
Department  of  Physics 
Los  Angeles  CA  90024>1547 


AFOSR-91-0006 


AEOSR- 


8.  PERFORMING  ORGANIZATION 
REPORT  number 

9  3  0^92 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 

AFOSR/NE 

110  Duncan  Avenue 
Bolling  AFB  DC  20333-0001 


11.  SUPPLEMENTARY  NOTES 

IZa.  DISTRIBUTION /AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

UNLIMITED 

This  docuzaent  has  been  apptoved 

for  public  release  and  sole;  itt 

13.  ABSTRACT  (Maximum  200  words) 

This  is  the  Final  Technical  Report  on  work  performed  under  the  si^port  of 
the  Air  Force  Office  of  Scientit  Research  under  Grant  AFOSR  91-0006  for 
the  peridd  October  1,  1992  to  Septeidaer  SO,  1993.  The  objective  of  this 
work  was  to  ceurry  out  basic  physics  research  on  the  generation  of  coherent, 
tunable  radiation  frcai  relativistic  electron  beams  and  exploring  meauis 
of  improving  performance. 


s 


DTIC 

ELECTE 
DEC  2  7  1993 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev  2-89) 

PFCKrib«d  by  ANSI  Sid  Z3V-18 
298-102 


Computer  Simulations  of  Radiation  Generation  From 
Relativistic  Electron  Beams 

AFOSR  91-0006  Final  Technical  Report 


Dr.  Anthony  T.  Lin 
UCLA  Physics  Department 

DTK?  OftlALmr  rwaPECTED  B 


University  of  California.  Los  Angeles 
Department  of  Physics 
Los  Angeles,  CA  90024-1547 


COMPUTER  SIMULATIONS  OF  RADIATION  GENERATION  FROM 
RELATIVISTIC  ELECTRON  BEAMS 

AFOSR  91-0006 


FINAL  TECHNICAL  REPORT 

TO 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 


Dr.  Anthony  T.  Lin 


Department  of  Physics 
University  of  California,  Los  Angeles 
Los  Angeles,  CA  90024-1547 


October  1, 1992,  to  September  30, 1993 


TABLE  OF  CONTENTS 


Page 


I.  Introduction . 3 

I I .  Summary  of  Work  Accomplished . 3 

A.  Peniotron  Forward  Wave  Oscillators . .4 

B.  A  Harmonic  FEL  Based  on  The  Second  Harmonic 

Gyroresonance . .9 

C.  Dependence  of  Efficiency  on  Magnetic  Field  in  Gyro 

Backward  Wave  Oscillators . 9 

D.  Phase  Locking  of  Gyro  Backward  Wave  Oscillators  .  .  .10 

III.  AFOSR  Supported  Publications  (October  1 , 1 992- 

September  30, 1993) . 13 

APPENDICES . 14 


2 


I.  introduction 


This  is  the  Final  Technical  Report  on  work  performed  under  the  suf^miI  of  the 
Air  Force  Office  of  Scientific  Research  under  Grant  AFOSR  91-0006  for  the  period 
October  1,  1992  to  September  30.  1993.  The  objective  of  this  work  was  to  carry  out 
basic  physics  research  on  the  generation  of  coherent,  tunable  radiation  from 
relativistic  electron  beams  and  exploring  means  of  improving  performances. 

II.  Summary  of  Work  Accomplished 

This  work  covered  computer  simulations  and  theoretical  studies  of  radiation 
generation  in  free  electron  microwave  devices.  We  have  carried  out  computer 
simulations  and  theoretical  analyses  to  address  various  issues  concerned  with  free 
eiectron  devices  which  includes  free  electron  lasers,  gyro  backward  wave  oscillators, 
and  peniotron  fonward  wave  oscillators. 

We  have  conceived  a  peniotron  forward  wave  oscillator  which  relies  on  the 
pulse  spreading  of  the  forward  propagating  absolute  instability  in  peniotron  interaction 
to  provide  the  feedback  mechanism.  Aside  from  the  usual  advantages  (high  efficiency 
at  high  harmonic  operation)  of  peniotron  interaction,  the  proposed  scheme  requires 
reiatively  short  interaction  length  and  can  couple  the  wave  power  out  from  the  beam 
exit  side  which  render  building  a  peniotron  oscillator  feasible  and  desirable. 

We  have  devised  a  new  concept  for  harmonic  free  electron  lasers  which  is 
based  on  our  early  theoretical  prediction  that  the  axial  guide  field  can  be  explored  as  a 
harmonic  selection  mechanism.  Computer  simulations  were  carried  out  which 
demonstrated  that  a  second  harmonic  FEL  based  on  the  gyroresonance  effect  can 
attain  gain  and  efficiency  comparable  to  the  fundamental  FEL. 

In  order  to  interpretate  the  recent  experimental  results  which  show  that  gyro 
backward  wave  oscillators  operated  close  to  the  waveguide  cutoff  can  achieve  much 
higher  efficiency,  computer  simulations  were  performed.  The  results  demonstrate  that 
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one  of  the  plausible  mechanisms  is  the  resulting  initial  large  frequency  mismatch 
(detuning)  which  places  the  electron  bunch  in  an  azimuthal  phase  close  to  where  the 
maximum  electron-wave  coupling  occurs. 

To  control  the  output  phase  of  an  osdliator.  an  external  signal  close  to  the  free- 
running  frequency  is  often  injected  to  phase  lock  the  oscillator.  In  carrying  out 
computer  simulations  to  illustrate  the  phase  locking  phenomenon,  an  unexpected 
effect  that  the  injected  signal  also  enhance  the  output  efficiency  is  observed. 

Below  we  summarize  the  scope  and  main  findings  of  these  investigations. 
More  detailed  information  are  contained  in  the  appended  reports. 

A.  EfinifltfQn  Forward  Wave  Oscillators  (AopfiDdLs.1) 

Microwave  source  based  on  peniotron  mechanism  was  conceived 
almost  three  decades  ago\  However,  up  to  now  there  is  still  no  convincing  evidence 
to  assure  that  high  efficiency  interaction  has  been  achieved  in  any  of  the  experiments 
reported.  From  the  past  experience  a  consensus  has  been  reached  that  the  essential 
elements  for  efficient  peniotron  operation  are  a  vane  waveguide  and  an  axis- 
encircling  beam. 

An  unavoidable  feature  of  a  vaned  waveguide  is  the  small  frequency  separation 
between  the  modes  of  larger  azimuthal  mode  number  n  which  renders  the  mode 
competition  in  peniotron  device  an  important  issue  to  be  addressed.  Using  the 
parameters  of  Table  I,  the  dispersion  curves  of  the  relevent  waveguide  modes  are 
shown  in  Fig.  1 .  The  first  three  cyclotron  harmonic  (s  =  1 ,  2,  3)  lines  are  also  plotted  in 
the  same  figure.  There  are  three  potential  oscillation  sources  which  are  marked  by 
points  1 ,  2  and  3.  They  all  have  their  own  starting  oscillation  lengths.  Points  2  and  3 
are  the  conventional  gyrotron  backward  wave  modes  with  relatively  low  efficiency. 
Point  1  is  the  peniotron  mode. 
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In  carrying  out  the  stability  analysis  of  a  peniotron  amplifier,  we  realized  that  the 
intrinsic  feedback  mechanism  in  an  absolute  instability  even  though  is  deleterious  to 
an  amplifier  can  be  utilized  to  form  oscillator  without  a  cavity  structure.  An  unstable 
system  will,  in  general,  be  unstable  for  a  range  of  kz  values  and  an  initial  disturbance 
will  excite  a  spectrum  of  kz's.  In  a  weakly  unstable  system  with  beam  current  below 
the  threshold  of  the  absolute  instability  (point  1).  the  instability  is  convective.  The 
disturbance  propagates  in  the  forward  direction  and  the  pulse  will  spread  to  each  side 
of  the  peak.  As  the  beam  current  is  increased  to  such  an  extent  so  that  the  growth  rate 
is  strong  enough  to  cause  the  signal  to  grow  at  the  very  point  of  the  initial  disturbance 
by  pulse  spreading,  the  instability  becomes  absolute.  As  is  known  that  an  absolute 
instability  grows  from  the  initial  fluctuation  of  the  system  eigenmode  and  is  eventually 
saturated  by  some  nonlinear  mechanism.  Therefore  by  properly  choosing  the 
parameters  one  should  be  able  to  turn  it  into  an  useful  microwave  source.  The  time 
evolution  of  the  wave  power  (TE31  mode)  from  the  peniotron  interaction  at  z  **  0  and  L 
(L  =  5  cm)  is  plotted  in  Fig.  2a.  The  results  clearly  illustrate  that  initial  perturbations 
grow  in  time  at  every  point  in  the  interaction  space  with  the  same  growth  rate  which  is 
the  signature  of  an  absolute  instability.  The  direction  of  power  flow  determined  by 
performing  local  Poynting  flux  calculation  x  is  along  the  beam  propagating 
direction  and  the  output  power  attained  more  than  50%  of  the  beam  input  power.  The 
electron  perpendicular  velocity  in  the  entire  interaction  space  at  ooct  =  6000  is 
displayed  in  Fig.  2b.  Simulation  results  clearly  reveal  that  as  a  result  of  the  peniotron 
interaction  all  electrons  irrespective  of  their  initial  phase  give  up  their  energies  to  the 
wave. 
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TABLE  I 

Second  Harmonic  Vaned  Peniotron  Oscillator  Design 


Beam  Voltage 

70  kV 

Beam  Current 

3.5  A 

Magnetic  Field 

6.5  kG 

Mode 

n  mode 

Number  of  Vanes 

6 

a  »  vj.  /  V|( 

1.2 

Bq/  Bg 

1 

0]c/2jc 

33.65  GHz 

Inner  Circuit  Radius  (a) 

0.16  cm 

Outer  Circuit  Radius  (b) 

0.304  cm 

6 


k^c 


Figure  1 :  Dispersion  cuves  of  a  six  vane  slotted  waveguide;  potential 
self  oscillation  modes  are  marked  by  1.  2  and  3. 
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B.  A  Harmonic  FEL  Based  on  The  Second  Harmonic  Gyroresonance 


(Appendix  2) 

In  recent  years,  there  is  interest  in  the  generation  of  harmonics  in  free  electron 
lasers  in  an  effort  to  achieve  ever  shorter  wavelengths  without  the  necessity  of  ever 
higher  electron  beam  energies  or  shorter  wiggler  periods.  The  conventional  approach 
for  generating  harmonic  emission  requires  a  strong  wiggler  field  to  produce  a  periodic 
longitudinal  motion  of  electrons  which  produces  gain  in  higher  harmonics.  In  this 
scheme,  a  substantial  part  of  beam  kinetic  energy  is  spent  in  exciting  the  fundamental 
emission  and  thus  the  harmonic  radiation  has  often  been  limited  by  problems  of 
insufficient  gain  and  low  efficiency.  Therefore,  harmonic  FEL  operation  generally 
requires  a  mechanism  that  selectively  enhances  the  competitiveness  of  a  desired 
harmonic.  We  propose  a  different  concept  based  on  our  early  theoretical  prediction^ 
that  the  axial  guide  field  can  be  explored  as  a  harmonic  selection  mechanism.  With 
the  guide  field  tuned  near  to  the  nth  wiggler  harmonic,  the  fundamental  harmonic 
component  of  the  quiver  velocity  will  be  reduced  by  a  factor  approximately  equal  to  n, 
whereas  its  nth-harmonic  component  is  resonantly  enhanced.  The  appropriate  beam 
source  for  such  an  application  is  preferably  annular  in  shape  with  a  resonably  large 
radius.  So,  upon  resonant  action  in  the  wiggler  region,  the  electrons  will  acquire  a 
substantial  harmonic  quiver  velocity  with  a  minimum  cross-sectional  velocity  spread. 
We  have  carried  out  computer  simulations  and  demonstrated  the  feasibility  of  a 
second  harmonic  FEL  based  on  the  gyrosonance  effect  with  attractive  gain  and 
efficiency.  Preliminary  results  also  indicate  a  susceptibility  to  interference  from  higher 
harmonic  modes. 

C.  Dependence  of  Efficiency  on  Magnetic  Field  in  Gvro  Backward 
Wave  Oscillators  (Appendix  3) 
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One  of  the  advantages  of  utilizing  gyro  BWO  as  millimeter  wave  source  is  its 
tunability  through  varying  the  magnetic  field.  Therefore  it  is  important  to  investigate 
how  the  variation  of  magnetic  field  influences  the  output  efficiency.  As  is  expected 
simulation  results  show  that  the  output  frequency  scales  almost  linearly  with  the 
magnetic  field  in  most  of  the  region  except  very  near  the  cutoff  while  the  output 
efficiency  exhibits  a  significant  jump  (about  50%)  when  the  oscillation  frequency 
approaches  the  waveguide  cutoff  frequency.  As  a  result  a  roughly  constant  output 
power  gyro  BWO  can  be  operated  in  the  following  two  regimes:  high  efficiency  with 
about  8%  tunability  or  medium  efficiency  with  about  40%  tunability. 

In  tjie  case  of  coo  =  1  043  coc  the  output  effieciency  is  about  20%.  This  is  in  good 
agreement  with  the  result  of  a  recent  experiment^.  In  order  to  illustrate  how  the 
electron  azimuthal  bunching  process  differs  when  the  magnetic  field  is  varied,  the 
phase  trajectories  of  forty  test  electrons  initially  (z  =  0)  uniformly  distributed  between  0 
and  2n  are  followed.  Simulation  results  show  that  one  of  the  plausible  mechanisms 
for  this  high  efficiency  operation  is  the  resulting  substantially  larger  frequency 
mismatch  (detuning)  which  places  the  electron  bunch  close  to  8  =  ^  where  the 

maximum  electron-wave  coupling  takes  place.  At  the  same  time  the  group  velocity  of 
the  wave  in  this  case  is  also  slower  which  increases  the  electron-wave  interaction 
time.  These  have  all  contributed  to  the  enhancement  of  output  efficiency. 

D.  Phase  Locking  of  Gvro  Backward  Wave  Oscillators 

Injecting  a  weak  signal  into  a  more  powerful  free-running  gyrotron  backward 
wave  oscillator  can  produce  interesting  and  useful  effects,  such  as  controlling  the 
phase  of  the  output  wave.  A  series  of  simulations  have  been  carried  out  to  study  the 
phase  locking  phenomenon  by  using  the  configuration  of  the  recent  experimental  set¬ 
ups.  The  scheme  for  diagnosing  whether  an  oscillator  has  been  phase  locked  or  not  is 
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by  mixing  the  injected  signal  (coj)  with  the  output  signal  (cdo)  and  taking  the  time 
averaging  over  one  injected  signal  period. 


I  <  Eo(t)  sinmjt  >  1 
1  <  Eo(t)  COSCDjt  >  J 


»Acoot-(j)i  .  (1) 

where  Acoo  s  Oj  -  cdq  and  <|>i  is  the  phase  difference.  The  ouutput  wave  is  phase  locked 
by  the  input  signal  if  <|>(t)  asymtotically  approaches  a  constant.  Figure  3  shows  the 
simulation  results  using  Acoo  =1  X  10’^  coo  and  P/Po  =  -15dB.  The  external  signal  was 
injected  after  the  output  power  reached  steady-state  (coct  «  4000).  Figure  3b  clearly 
illustrates  that  <t>(t)  asymptotically  approaches  a  constant  which  is  the  signature  of 
phase  locking.  The  conventional  theory  predicts  that  the  output  power  can  only  be 
increased  by  the  amount  of  the  injected  signal  power.  Figure  3a  reveals  that  in  the 
presence  of  the  injected  signal  the  output  power  is  enhanced  by  almost  50%.  This 
unexpected  effect  requires  more  study.  In  contrast  to  the  classical  oscillator,  the 
quality  factor  of  a  gyro  backward  wave  oscillator  is  undefined  which  could  complicate 
the  theoretical  analysis. 


4000 


6000  8000 

Wet 


Figure  3:  Gyro  backward  wave  oscillator  with  an  external  signal: 
(a)  output  pulse  and  (b)  showing  injection  locking. 
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Abstract 

The  pulse  spreading  of  the  forward  propagating  absolute  instability  in  gyro- 
peniotron  interaction  can  be  utilized  to  provide  the  feedback  mechanism  for 
establishing  forward  wave  oscillators.  By  properly  choosing  the  interaction  length  and 
beam  current  high  efficiency  stable  output  at  high  harmonic  operation  can  be  attained. 
This  is  demonstrated  by  performing  computer  simulations  which  also  reveal  that  the 
peniotron  interaction  is  able  to  suppress  the  gyrotron  interaction  and  reaches  the 
optimal  efficiency  without  the  influence  of  the  gyrotron  mode. 
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L  Introduction 

Microwave  source  derived  from  peniotron  mechanism  was  conceived  almost 
three  decades  ago^  Since  then  many  investigators  have  carried  out  theoretical 
calculations^'^  and  attempted  to  employ  various  beam  and  electrodynamic 
configurations*'^  to  illustrate  the  aspect  of  high  efficiency  peniotron  harmonic 
interaction.  Up  to  now  there  is  still  no  convincing  evidence  to  assure  that  high 
efficiency  interaction  has  been  achieved  in  any  of  the  experiments.  However  from  the 
past  experience  a  consensus  has  been  reached  that  the  essential  elements  for 
efficient  peniotron  operation  are  a  electrodynamic  configuration  which  is  capable  of 
providing  a  positive  radial  gradient  in  the  microwave  field  distribution  and  an  axis- 
encircling  beam  which  initially  can  be  placed  at  a  high  beam-wave  coupling  position. 
Recently,  based  on  these  important  informations,  the  parameters  of  a  peniotron 
amplifier*  were  proposed  and  a  high  quality  factor  (Q  z  2000)  peniotron  oscillator  was 
built^.  The  former  has  to  avoid  various  self  osdilations  which  could  degrade  or  totally 
prohibit  the  amplification  process.  This  constraint  forces  one  to  contemplate  a  multi¬ 
section  device  with  the  length  of  each  section  below  the  starting  oscillation  length  of  all 
the  potential  oscillations.  The  high  quality  factor  of  the  latter  tends  to  require  a  well 
matched  external  load  to  couple  out  the  major  portion  of  the  stored  energy.  It  is  also 
suspected  that  the  experimental  peniotron  output^  may  be  contaminated  by  the  nearby 
gyrotron  oscillation.  In  this  paper  we  propose  an  alternative  peniotron  oscillator 
configuration  without  external  reflection  which  alleviates  the  difficulty  of  finding  a  well 
matched  load  circuit.  The  feedback  mechanism  is  provided  by  the  pulse  spreading  of 
the  forward  propagating  absolute  instability  which  can  be  chosen  to  be  the  dominant 
interaction  in  the  system.  The  proposed  configuration  may  be  considered  as  the  first 
step  to  demonstrate  that  peniotron  interaction  is  capable  of  attaining  high  efficiency,  in 
investigating  the  peniotron  and  gyrotron  competition,  simulation  results  reveal  that  the 
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peniotron  interaction  always  suppresses  the  gyrotron  interaction  if  the  beam  guiding 
center  spread  is  moderate  (<  30  %).  This  understanding  should  be  helpful  in 
interpretating  the  experimental  result  of  conventional  peniotron  oscillators^. 
Simulation  results  are  also  capable  of  determining  the  starting  oscillation  length  of 
absolute  instability  which  is  an  important  information  for  a  peniotron  amplifier. 

II.  Numerical  Model 

Because  of  its  boundary  conditions,  the  harmonic  field  in  an  azimuthally 
corrugated  interaction  structure  originally  introduced  for  high  harmonic  gyrotron 
operations^' becomes  stronger  when  it  gets  closer  to  the  waveguide  wall  and  is 
therefore  also  suitable  for  the  peniotron  interaction.  The  cross-sectional  view  of  a  six 
vane  slotted  waveguide  with  an  axis-encircling  beam  is  displayed  in  Fig.  1a.  In  the 
interaction  region,  the  electric  and  magnetic  fields  of  a  transverse  electric  field 
mode^^'^^  may  be  written  as 


En  =  Ejn(Z,t)  ( 2  X  Vj^Cn)  (1 ) 

Bn  =  BTn(2,t)  Vj^Cn  +  2  Bi^pk^Cn  (2) 

where  Cn  is  the  local  wave  function 

oo 

C„.-  ^^Jr{kxr)el™  .  (3) 

m=-oo 


Here  r  =  n  +  mN,  n  =  0,1 . N/2,  m  is  any  integer,  N  is  the  number  of  vanes,  and 
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Ar  »  is  the  enhancement  factor  of  the  slotted  structure,  kj.  is  the  cutoff 

Ji  rSo 

wavenumber  which  can  be  determined  from  the  circuit  dispersion  relation 


Jr(Kia)  /  sinrep  \2 
Jr'(kia) '  rep  ' 


n  Zp(kj.a) 

N0pZi(kxa)  ’ 


(4) 


where  Zr  is  a  combination  of  Bessel  function  (Jr)  and  Neumann  function  (Np)  of 
order  r 


Zr(x)=Jr(x)-i;j;^Nr(x)  .  (5) 

and  do=2n  /  4N.  Here  a  and  b  denote  respectively  the  inner  and  outer  circuit  radius. 
For  each  value  of  n,  there  are  infinite  number  of  solutions  from  Eq.(4).  The  first  root  is 
the  most  important  one  and  the  associated  mode  is  usually  designated  by  TEni-  The 
mode  number  n  is  also  the  number  of  times  the  RF  field  pattern  repeats  in  one  rotation 
around  the  waveguide  axis.  There  are  many  possible  modes  but  the  two  most 
interesting  cases  are  the  n  mode  where  adjacent  slots  are  out  of  phase  by  n,  and  the 
2n  mode  where  the  phase  in  each  slot  is  identical.  In  the  most  rigorous  description  an 
infinite  superposition  of  different  Fs  is  required  to  represent  each  n  mode.  However  in 
most  of  the  cases  an  adequate  approximation  includes  just  two  partical  waves  -  those 
with  r  =  n  and  (n  -  N).  The  physical  meaning  of  this  is  that  an  electromagnetic  wave 
rotating  in  a  periodic  structure  must  interfer  with  other  waves  to  reproduce  the  static 
periodic  pattern  of  the  structure. 

The  modes  supported  by  the  slotted  structure  can  interact  with  an  axis- 
encircling  electron  beam  through  the  relativistic  gyrotron  and  non-relativistic  peniotron 
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mechanisms.  In  a  gyrotron  interaction  it  is  necessary  that  an  electron  stays 
approximately  in  the  same  phase  with  respect  to  the  wave  as  it  traverses  through  the 
interaction  region  so  that  the  cumulative  electron-wave  interaction  produces  a 
substantial  bunching  through  relativistic  mass  effect.  As  a  consequence,  an  initially 
randomly  phased  electron  beam  would  end  up  with  some  of  its  electrons  giving  up 
energy  to  and  the  rest  of  them  absorbing  energy  from  the  same  wave.  This  resonant 
condition  can  be  expressed  as 


0)0  =  kjVz  +  s 


Oc 


(6) 


where  o)o  and  kj  are  respectively  the  wave  frequency  and  axial  wavenumber,  vj  and  y 
are  the  electron  axial  velocity  and  relativistic  factor.  Qc  is  the  non-relativistic  cyclotron 
frequency  and  s  is  the  harmonic  number.  For  efficient  gyrotron  interaction  n  must  be 
equal  to  s.  On  the  other  hand  the  gain  mechanism  underlines  a  peniotron  interaction 
is  due  to  the  electron  guiding  center  drift  so  that  during  a  cyclotron  orbit  an  electron  on 
the  average  is  in  a  stronger  RF  field  region  during  the  decelerating  phase  and  in  a 
weaker  RF  field  region  during  the  accelerating  phase.  This  is  due  to  the  fact  that  an 
accelerated  electron  increases  its  radius  of  gyration  and  the  requirement  n  =  s  +  1 
which  ensures  that  in  our  example  (Fig.  la)  half  of  the  cyclotron  orbit  is  accelerating 
while  the  other  half  Is  decelerating.  As  a  result  all  electrons  regardless  of  their  initial 
phases  relative  to  the  wave  convert  their  kinetic  energy  into  wave  energy  and  makes 
peniotron  devices  inherently  high-efficiency. 

An  unavoidable  feature  of  a  vaned  waveguide  is  the  small  frequency  separation 
between  the  modes  of  larger  n  which  renders  the  mode  competition  in  peniotron 
device  an  important  issue  to  be  addressed.  Using  the  parameters  of  Table  I,  the 
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dispersion  curves  of  the  relevent  waveguide  modes  are  shown  in  Fig.  1b.  The  first 
three  cyclotron  harmonic  (s  s  i ,  2,  3)  lines  are  also  plotted  in  the  same  figure.  There 
are  three  potential  oscillation  sources  which  are  marked  by  points  1 , 2  and  3.  They  all 
have  their  own  starting  oscillation  lengths.  Points  2  and  3  are  the  conventional 
gyrotron  backward  wave  modes  with  relatively  low  efficiency.  Point  1  is  the  desired 
peniotron  mode  which  is  capable  of  producing  high  efficiency  output  as  well  as  large 
temporal  growth  rate  due  to  its  grazing  interaction. 

III.  Simulation  Results 

The  intrinsic  feedback  mechanism  in  an  absolute  instability  even  though  is 
deleterious  to  an  amplifier  can  be  utilized  to  form  an  oscillator  without  a  cavity 
structure.  Gyrotron  backward  wave  oscillators^ which  rely  on  a  backward 
propagating  electromagnetic  wave  to  feedback  the  energy  generated  by  the  forward 
propagating  beam  have  been  built  and  performed  rather  well.  As  far  as  the  absolute 
unstable  mode  at  point  1  (Fig.  1b)  is  concerned,  experiments  carried  out  so  far  were 
designed  to  avoid  the  excitation  of  this  mode  by  using  parameters  below  its  starting 
oscillation  condition. 

The  theoretical  investigations  of  the  absolute  instability  for  either  gyrotrons^®*’® 
or  peniotrons^^  are  all  based  on  the  pinch-point  theory  which  was  first  derived  in 
plasma  physics  for  an  infinite  homogeneous  medium^®.  However,  an  absolute 
instability  in  a  limited  spatial  extent  may  be  stable.  In  order  to  illustrate  the  effects  of 
system  length  on  the  competition  among  various  absolute  unstable  modes,  we  have 
first  carried  out  single  transverse  mode  computer  simulations  to  evaiute  their  temporal 
growth  rate  versus  interaction  length  and  results  are  shown  in  Fig.  2a.  The  starting 
oscillation  length  can  be  determined  by  extending  the  growth  rate  curve  to  intersect 
the  z  axis.  The  time  evolution  of  the  wave  power  (TE31  mode)  from  the  peniotron 
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interaction  at  z  «  0  and  L  for  the  L  «  5  cm  case  is  plotted  in  Fig.  2b.  The  results  clearly 
illustrate  that  initial  perturbations  grow  in  time  at  every  point  in  the  interaction  space 
with  the  same  growth  rate  which  is  the  signature  of  an  absolute  instability.  The 
direction  of  power  flow  determined  by  performing  local  Poynting  flux  calculation 
Et  X  Bt  is  along  the  beam  propagating  direction  and  the  output  power  attained  more 
than  50%  of  the  beam  input  power,  it  was  also  observed  in  simulations  that  a  broad 
spectrum  of  k's  was  excited  and  they  all  propagate  in  the  forward  direction.  In  the 
range  of  L  ^  10  cm,  there  is  no  trace  of  exciting  the  backward  wave  of  point  3  (Fig.  1b). 
Using  the  same  system  length  and  retaining  only  the  TE21  mode,  the  time  evolution  of 
the  wave  power  at  both  ends  is  shown  in  Fig.  2c.  In  this  case  the  instability  arises  from 
the  gyrotron  backward  wave  interaction.  The  efficiency  has  reached  only  about  10% 
and  the  wave  power  flow  is  in  the  backward  direction.  Notice  that  the  wave  amplitude 
at  z  s  L  is  not  negligible  in  comparison  with  that  of  at  z  =  0  because  point  2  (Fig.  1b)  is 
so  close  to  the  waveguide  cutoff.  In  order  to  demonstrate  the  difference  between  the 
peniotron  and  the  gyrotron  interactions,  their  respective  electron  perpendicular 
velocity  in  the  entire  interaction  space  at  coct  s  6000  is  displayed  in  Fig.  3.  Simulation 
results  clearly  reveal  that  as  a  result  of  the  peniotron  interaction  all  electrons 
irrespective  of  their  initial  phase  give  up  their  energies  to  the  wave  while  the  gyrotron 
interaction  results  in  that  only  part  of  electron  beam  loses  energy  and  the  remaining 
part  gains  energy. 

An  unstable  system  will,  in  general,  have  a  range  of  kz  values  with  positive 
(Oi(kz)  and  an  initial  disturbance  will  excite  a  spectrum  of  kz's.  In  a  weakly  unstable 
system  with  beam  current  below  the  threshold  for  oscillation,  the  instability  is 
convective.  The  disturbance  like  point  1  in  Fig.  1b  propagates  in  the  forward  direction 
and  the  pulse  will  spread  to  each  side  of  the  peak.  As  the  beam  current  is  increased  to 
such  an  extent  so  that  the  growth  rate  is  strong  enough  to  cause  the  signal  to  grow  at 
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the  very  point  of  the  initial  disturbance  by  pulse  spreading,  the  instability  becomes 
absolute. 

The  threshold  of  the  absolute  instability  in  an  infinite  system  can  be  evaluated 
by  using  the  expression  derived  in  Ref.  [19]  which  describes  the  spatial  and  temporal 
evolution  of  the  response  electric  field  to  a  localized  initial  disturbance 

-  (2ni<»”(ko)tl'« 

where  o)=o)(k)  is  the  dispersion  relation  for  the  peniotron  interaction^ ^,Vg  = 
(Or'(ko)  where  ko  is  such  that  coi’(ko)  *  0  and  o)i"(ko)  <  0  and  the  prime  is  the 
derivative  with  respect  to  kz.  Equation  (7)  reveals  that  the  evolving  pulse 
peaks  at  z^Vgt,  the  pulse  shape  is  growing  and  spreading  in  time  if  a)j"(ko) 
is  negative,  and  its  mean-square  spatial  spread  is  the  reciprocal  of 
lm[^-7"— The  threshold  can  be  determined  from  the  fact  that  E(z»0,t) 

must  increase  exponentially  with  time  or 

lm[o)(ko)]  >  Im  [  g  •  (®) 

This  expression  only  gives  an  approximate  threshold  for  self  oscillation  in  an  infinite 
system.  By  substituting  the  parameters  used  in  simulations  into  the  dispersion  relation 
given  in  Ref.  [17],  we  obtained  respectively  0.0157  (Dc  for  the  left  hand  side  and  0.014 
(Dc  for  the  right  hand  side  of  Eq.  (8).  Clearly,  the  inequality  in  Eq.  (8)  is  satisfied.  The 
evidence  of  spatial  pulse  spreading  at  early  time  can  also  be  infered  indirectly  from 
Fig.  4a  which  displays  the  time  evolution  of  unstable  pulse  in  wave  number  space 
before  the  stored  energy  begins  to  exponentiate  in  time.  Notice  that  the  pulse  width  in 
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real  and  kj  space  is  inversely  proportional  to  each  other.  The  width  of  the  unstable  Kz 
spectrum  remains  relatively  constant  during  the  wave  growing  period  (Fig.  4b).  As  is 
known  that  an  absolute  instability  grows  from  the  initial  fluctuation  of  the  system 
eigenmodes  and  is  eventually  saturated  by  some  nonlinear  mechanism.  Therefore  by 
properly  choosing  the  parameters  one  should  be  able  to  turn  it  into  an  useful 
microwave  source.  Based  on  our  prelimary  simulation  results  which  will  be  presented 
later,  the  proposed  mechanism  possesses  certain  unique  characteristics  which  make  it 
a  very  attractive  and  viable  peniotron  microwave  source. 

To  assure  that  the  oscillation  from  the  TE21  mode  would  not  deteriorate  the 
performance  of  the  desired  peniotron  operation,  multimode  simulations  were  carried 
out.  Figure  5  shows  the  time  evolution  of  the  stored  energy  in  the  interaction  space  for 
three  different  cases.  With  L  =  5  cm  (Fig.  5a)  the  initial  growth  of  the  TE21  mode  comes 
from  the  convective  unstable  mode  with  positive  kz  and  the  growth  rate  of  its  absolute 
unstable  mode  is  so  small  that  it  can  be  easily  suppressed  by  the  growth  of  the  TE31 
mode.  Lengthening  the  system  length  to  8  cm  tends  to  enhance  the  gyrotron 
interaction  more  than  to  the  peniotron  interaction,  in  this  case  the  two  interactions 
have  about  the  same  growth  rate  (Fig.  5b).  The  peniotron  mode  eventually  is  able  to 
suppress  the  gyrotron  mode.  Using  L  =  8  cm  and  I  =  1.5  A  (Fig.  5c)  the  temporal 
growth  rate  of  the  gyrotron  interaction  now  surpasses  that  of  the  peniotron  interaction. 
However  the  peniotron  interaction  is  still  able  to  suppress  the  gyrotron  interaction.  The 
respective  time  evolution  of  the  peniotron  output  power  of  the  three  cases  described  in 
Fig.  5  is  shown  in  Fig.  6.  Both  single  mode  and  two  modes  simulations  are  presented. 
The  results  indicate  that  the  presence  of  the  gyrotron  interaction  when  it  reaches  large 
amplitude  may  affect  the  temporal  growth  rate  of  the  peniotron  interaction  which 
requires  no  phase  bunching  but  not  its  final  efficiency  for  all  three  cases.  This  is 
because  that  to  a  peniotron  interaction  a  saturated  gyrotron  mode  induces  effective 


beam  guiding  center  and  perpendicular  velocity  spreads  which  would  reduce  the 
temporal  growth  rate.  Since  at  the  end  the  gyrotron  mode  is  completely  suppressed 
and  returns  to  its  noise  level,  the  peniotron  mode  should  be  able  to  attain  the  efficiency 
without  the  influence  of  the  gyrotron  mode.  The  effects  of  guiding  center  spread  and 
beam  current  on  the  performance  of  the  oscillator  (1=5  cm)  are  respectively  displayed 
in  Fig.  7a  and  Fig.  7b.  The  results  show  that  a  20%  guiding  center  spread  would 
degrade  the  output  efficiency  to  40%  and  I  =  3.5  A  is  the  optimal  beam  current  for 
maximizing  the  effieciency.  The  optimal  system  length  (Lo)  which  maximizes  the 
output  efficiency  is  very  close  to  5  cm  (I  =  3.5  A).  With  L  >  Lo.  most  of  the  electrons  (in 
the  z  >  Lo  region)  due  to  changes  in  their  energy  and  guiding  center  location  resulting 
from  electron-wave  interactions  are  no  longer  axis  encircling  and  eventually  becomes 
resonance  with  the  wave.  As  a  result  of  detuning,  electrons  on  the  average  start  to 
absorb  energy  from  the  wave  and  terminates  the  wave  growth  process. 

IV.  Summary 

In  summary,  we  have  proposed  a  scheme  for  establishing  peniotron  forward 
wave  oscillator  without  external  reflection.  The  feedback  mechanism  relies  on  the 
pulse  spreading  of  the  absolute  instability.  Simulation  results  also  reveal  that  the 
peniotron  interaction  is  capable  of  suppresing  the  gyrotron  interaction.  This 
suppression  phenomenon  should  also  occur  in  a  conventional  peniotron  oscillator 
which  requires  less  beam  current  to  start  the  interaction. 
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TABLE  I 

Second  Harmonic  Vaned  Peniotron  Oscillator  Design 


Beam  Voltage 

70  kV 

Beam  Current 

3.5  A 

Magnetic  Field 

6.5  kG 

Mode 

n  mode 

Number  of  Vanes 

6 

a  =  vi  /  V|| 

1.2 

Bo/Bg 

1 

(i)c/27C 

33.65  GHz 

Inner  Circuit  Radius  (a) 

0.16  cm 

Outer  Circuit  Radius  (b) 

0.304  cm 
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Figure  Captions: 

Figure  1  Six  vane  siotted  waveguide  and  an  axis  encircling  beam:  (a)  cross-sectional 
view  and  (b)  dispersion  curves. 

Figure  2  Self  oscillations  in  finite  length  system:  (a)  temporal  growth  rate  versus 
system  length  (star:  TE31,  circle:  TE21)  and  the  time  evolution  of  wave 
power(Po)  normalized  to  beam  power  (Pb)  for  L  «  5  cm  (b)  peniotron 
interaction  (TE31 )  and  (c)  gyrotron  interaction  (TE21 ). 

Figure  3  Electron  perpendicular  velocity  versus  axial  distance:  (a)  peniotron 
interaction  and  (b)  gyrotron  interaction. 

Figure  4  Time  evolution  of  the  unstabie  kz  spectrum:  (a)  before  the  stored  energy 
begins  to  grow  and  (b)  during  the  growth. 

Figure  5  Time  evolution  of  stored  wave  energy  from  multi  mode  simulations,  (a)  L  *  5 
cm,  I  s3.5  A,  (b)  L  =  8  cm,  I  =  3.5  A,  and  (c)  L  =  8  cm,  I  *  1.5  A. 

Figure  6  Time  evolution  of  peniotron  output  wave  power  (solid  curve  is  from  single 
mode  simulation,  dashed  curve  is  from  multi  mode  simulation  :  (a)  L  s  5  cm, 
I  =3.5  A,  (b)  L  =  8  cm,  I  =  3.5  A,  and  (c)  L  =  8  cm,  I  =  1 .5  A. 

Figure  7  Effects  of  (a)  guiding  center  spread  and  (b)  beam  current  on  the 
performance  of  the  oscillator. 
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A  harmonic  FEL  based  on  the  second  harmonic  gyroresonance 
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The  solenoidal  guide  field  of  Raman  free  electron  lasers  can  be  tuned  to  resonantly  enhance  a  desired  harmonic  component  of 
the  electron  quiver  velocity.  The  feasibility  of  a  harmonic  FEL  based  on  this  effect  is  investigated.  Computer  simulations 
demonstrate  that  the  gain  and  efficiency  of  the  second  harmonic  generation  can  be  comparable  to  that  of  the  fundamental 
generation. 


1.  Introduction  2.  Harmonic  gyroresonance 

Off-axis  electron  motion  is  not  amenable  to  exact 
In  order  to  transport  the  high  current  beam  (a  few  analysis.  It  has  been  studied  numerically  [3-5]  or  by 

hundred  A)  through  the  wiggler  region  of  a  Raman  some  averaging  method  [5,6).  In  these  treatments,  the 

FEL,  an  axially  directed  guide  magnetic  field  is  often  harmonic  aspect  of  the  electron  motion  has  been  either 

employed  to  counter  the  defocusing  effect  of  the  obscure  J  or  averaged  out.  In  ref.  [7]  we  focused  on  the 

space-charge  field.  The  guide  field  exerts  a  strong  harmonic  content  of  the  election  motion  by  ignoring 

influence  on  the  electron  equilibrium  orbit  [1-7]  by  other  first  order  complexities  of  the  orbital  motion, 

which  the  gain,  the  operating  frequency,  and  the  out-  gj  ,(,g  betatron  oscillation  [5,6]  and  guiding  center 

put  efficiency  of  a  Raman  FEL  is  affected.  The  impor-  [34].  that  model  [7],  the  magnetic  forces  of  an 

tance  of  the  guide  field  to  the  performance  of  an  ^L  off-axis  electron  is  evaluated  on  an  ideal  helical  orbit, 

was  clearly  demonstrated  in  a  recent  FEL  amplifier  gj  projected  in  the  cross-sectional  plane  in  fig.  la. 

experiment  [8]  in  which  the  FEL  efficiency  was  in-  jbe  wiggler/guide  field  is  expressed  in  the  cylindrical 

creased  from  2%  to  27%  just  by  reversing  the  guide  coordinates  (r,  0,  z)  by  [2] 

field  direction.  Recent  multimode  computer  simulation 

results  [9]  suggest  a  possible  interpretation  of  the  dra-  jj  - -1-  2R*  cos(»  -  k^z)e, 

matic  contrast  in  efficiency:  the  normal  field  FEL 

provides  more  favorable  conditions  for  superradiant  /,(A:«r) 

amplification  of  local  noise  which  degrades  the  beam  sin(9  - 

quality  and  accounts  for  the  much  lower  gain  and  * 

efficiency.  -l-/,(fe„r)  sin(«  -  ,  (1) 

The  same  experiment  also  revealed  that  in  addition 

to  the  commonly  observed  dip  in  output  power  when  where  -  2it/A„,  I„  is  the  modified  Bessel  function 

the  wiggler  wavelength  (A.)  is  close  to  the  pitch  (Ap  =  of  the  nth  first  kind,  and  I’  is  its  derivative.  In  fig.  la 

of  the  electron  gyration  in  the  axial  guide  the  electron  phase  angle  is  chosen  so  that  in  the  limit 

field,  an  unexpected  dip  occurs  when  A,,  *  -Ap.  This  _»  q,  the  orbit  reduces  to  the  exact  steady-state  orbit, 

was  shown  [7]  to  be  the  effect  of  harmonic  gyroreso-  Defining  and  as  unit  vectors  perpendicular 

nance  of  the  off-axis  electrons.  Selective  enhancement  and  tangential  to  the  circle  of  fig.  la,  respectively,  we 

of  a  particular  harmonic  component  of  the  electron  niay  write  the  transverse  component  of  the  magnetic 

quiver  velocity  by  the  guide  field  suggests  the  interest-  as 

ing  possibility  of  harmonic  FEL  interactions.  In  this  j 

paper,  a  second  harmonic  FEL  based  on  the  harmonic  F ^  +  l2(,k,^r)  cos  26,]*^ 

gyroresonance  effect  is  demonstrated  through  com-  ^  r  \  ■  \ 

puter  simulations.  +^2(*w'')  26,e^|,  (2) 
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Table  1 

Paramelers  and  results  of  computer  simulations 


V 

(beam  energy) 

700  keV 

y 

(relativistic  factor) 

2.37 

P 

(t'/c) 

0.907 

I 

(beam  current) 

200A 

r. 

(initial  beam  radius) 

0.128  cm 

r* 

(waveguide  radius) 

0.255  cm 

<Vc/2ir 

(waveguide  cutoff  frequency) 

34.5  GHz 

A, 

(wiggler  period) 

3.18  cm 

Bo 

(guide  field) 

14.2  kG 

B, 

(wiggler  field) 

1  kG 

fo 

(input  frequency) 

60  GHz 

r 

(gain) 

1.2  dB/cm 

V 

(efficiency) 

6% 

kc 


Fig.  2.  Collective  interaction,  (a)  Dispersion  curves  of  the 
waveguide  and  beam  modes,  (b)  the  spatial  growth  profile  of 
the  injected  signal,  (c)  the  frequency  spectrum  at  the  output 
plane. 


this  cotiHguration,  only  the  second  harmonic  beam  line 
is  capable  of  interacting  with  the  waveguide  mode. 

As  is  shown  in  eq.  (9),  the  harmonic  content  of  the 
electron  quiver  velocity  is  of  the  order  of 
for  <  1.  In  order  to  conflrm  this  scaling,  the  orbits 
of  two  electrons  initially  located  at  r, » O.ZSr,^  and 
O.Sr^  were  followed  in  time  with  the  magnetic  field 
tuned  at  the  second  harmonic  resonance.  The  y-com- 
ponent  of  the  quiver  velocity  is  displayed  in  w-space  in 
figs,  lb  and  Ic.  The  electron  further  away  from  the  axis 
acquired  a  much  stronger  harmonic  component  while 
at  the  same  time  nearly  grazed  the  waveguide  wall. 
Therefore,  we  have  chosen  for  our  simulation  an  annu¬ 
lar  beam  with  «  O.Sr,..  Using  this  beam,  we  carried 
out  computer  simulations  of  an  FEL  amplifier  with  a 
system  length  of  eight  wiggler  wavelengths.  The  spatial 
growth  of  a  100  kW  injected  signal  at  (i>o“  I.73a>^  is 
shown  in  fig.  2b.  The  gain  and  efficiency  estimated 
from  simulation  results  are  respectively  1.2  dB/cm  and 
6%  which  are  comparable  to  that  of  the  fundamental 
FEL  amplifier.  The  spatial  profile  of  the  amplified 
signal  also  exhibits  a  regular  beating  structure.  Spectral 
analysis  carried  out  at  the  output  plane  (fig.  2c)  reveals 
that  in  addition  to  the  injected  signal  there  are  a  few 
higher  frequency  modes  (at  substantial  levels)  with 
frequency  separation  of  Aw  « nk^^i  g  from  the  main 
mode.  The  interference  of  these  m^es  with  the  in¬ 
jected  signal  results  in  the  observed  spatial  field  pat¬ 
tern  and  lower  efficiency  in  comparision  with  the  single 
mode  simulation. 

In  conclusion,  we  have  demonstrated  the  feasibility 
of  a  second  harmonic  FEL  based  on  the  gyroresonance 
effect  with  attractive  gain  and  efficiency  Preliminary 
results  also  indicate  a  susceptibility  to  interference 
from  higher  harmonic  modes. 
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A  gyro-backward  wave  oscillator  with  its  oscillation  frequency  close  to  the  waveguide  cutoff 
frequency  has  been  experimentally  demonstrated  by  other  investigators  to  attain  near  20% 
efficiency.  Computer  simulations  have  been  carried  out  to  elucidate  the  mechanisms  which  give 
rise  to  t^  high  efficiency  regime.  The  result  shows  that  one  of  the  plausible  mechanisms  is  the 
resulting  large  frequency  mismatch  (detuning)  which  places  the  electron  bunch  in  azimuthal 
phase  close  to  where  the  maximum  electron-wave  coupling  occurs. 


I.  INTRODUCTION 

The  gyrotron  is  a  mature  subject''^  and  has  been  rou¬ 
tinely  built  at  various  industrial  companies  for  the  use  of 
electron  cyclotron  resonant  heating  in  magnetic  fusion  de¬ 
vices.  In  order  to  achieve  high  efficiency,  it  often  requires 
the  magnetic  field  to  satisfy  grazing,  condition  which  limits 
its  tunability.  By  raising  the  sole  idal  magnetic  field,  the 
fundamental  cyclotron  line  will  eventually  intersect  the 
waveguide  mode  in  the  negative  k,  direction  of  &  to  (fre¬ 
quency)  versus  (axial  wave  number)  diagram  (Fig.  1). 
In  this  region  an  electron  beam  carries  its  energy  forward 
while  the  waveguide  mode  transports  the  generated  wave 
energy  backward.  This  provides  an  internal  feedback 
mechanism  and  renders  a  cavity  structure  unnecessary. 
Furthermore,  the  relaxation  of  the  grazing  condition 
makes  a  gyro-backward  wave  oscillator  (BWO)  tunable  by 
varying  either  the  beam  energy  or  magnetic  field.  Due  to  its 
tunability  and  simple  configuration,  gyro-BWO  has  re¬ 
cently  been  envisaged  as  a  viable  millimeter  wave  somce.^ 

In  the  past  there  were  a  few  theoretical  analyses  which 
determined  the  critical  oscillation  length  of  a  gyro-BWO.^’^ 
Lately,  the  theoretical  approach  had  been  improved  to  en¬ 
able  it  to  predict  the  eigenmode  temporal  growth  rate  of 
the  oscillation.^  An  up-tapered  magnetic  field  was 
observed^  to  increase  the  efficiency  from  around  10%  to 
30%  and  the  mechanism  responsible  for  the  enhancement 
was  recently  identified.'  Early  experimental  results^  were 
able  to  obtain  moderate  efficiency  but  disappointing  output 
power  which  was  limited  by  the  beam  current.  A  recent 
experiment‘d  employing  high-energy  electron  beams  was 
able  to  produce  high-power  (tens  of  megawatts)  but  low 
efficiency.  Very  recently,  it  was  demonstrated' '  experimen¬ 
tally  that  operating  close  to  the  waveguide  cutoff  fre¬ 
quency,  the  efficiency  of  a  gyro-BWO  can  reach  around 
20%  with  more  than  1(X)  kW  output  power.  In  this  paper 
we  shall,  through  computer  simulations,  offer  a  plausible 
explanation  for  the  efficiency  dependence  on  the  solenoidal 
magnetic  field  in  gyro-backward  wave  oscillators. 


The  efficiency  of  the  electron-wave  interaction  process  de¬ 
pends  on  the  phase  of  the  bunched  electrons  relative  to 
electromagnetic  wave.  To  investigate  how  the  change  of 
magnetic  field  affects  the  device  efficiency  we  shall  start 
with  a  brief  derivation  of  the  equation  governing  the 
electron-wave  energy  exchange.  Since  a  circularly  polar¬ 
ized  wave  is  a  good  approximation  of  the  TEn  mode, 
which  is  the  mode  we  are  interested  in,  we  assume  that  the 
electric  field  of  the  wave  can  be  written  as 


<0 

E=— lAUxcos^-f-ysin^),  (1) 

c 

where  c  is  the  speed  of  light,  A  is  the  vector  potential,  and 
is  the  phase  of  the  backward  propagating  wave. 
The  perpendicular  velocity  component  of  an  electron  mov- 
ing  in  a  uniform  magnetic  field  can  be  expressed  as 

sizVi  (.X  cos  i(>+y  sin  il>),  (2) 

where  ^  is  the  electron  gyrophase  angle.  Define 


(3) 

as  the  electron  phase  relative  to  the  wave.  The  addition  of 
ir/2  in  the  expression  for  0  is  purely  for  the  convenience  of 
simulation  diagnoses.  The  equation  describing  the  time  ev¬ 
olution  of  electron  energy  (y)  when  it  interacts  with  an 
electromagnetic  wave  is 


dt 


— evE 


me 


-sz—aaPx  sin  0, 


(4) 


where  m  is  the  electron  rest  mass,  fl=c|Al/mc^,  and 
/3^  =Vx  /c.  Equation  (4)  clearly  reveals  that  an  electron 
loses  energy  when  ir>0>O  and  gains  energy  when 
2ir>0> IT.  The  conversion  of  electron  energy  into  wave 
energy  is  most  efficient  when  0=ir/2. 

The  resonant  condition  between  a  backward  propagat¬ 
ing  electromagnetic  wave  (ci>o,A:o)  and  an  electron  in  an 
axial  magnetic  field  is 


II.  ELECTRON  PHASE  RELATIVE  TO 
ELECTROMAGNETIC  WAVE 

The  dominant  bunching  mechanism  in  electron  cyclo¬ 
tron  masers  (gyro  devices)  is  azimuthal  phase  bunching. 


where  fl,  is  the  nonrelativistic  electron  cyclotron  fre¬ 
quency.  The  electron  phase  versus  distance  is  therefore 
roughly  governed  by  the  following  equation: 
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FIG.  1.  The  tt)  vs  A,  diagram  for  a  gyro-backward  wave  oscillator  with 
K-50  kV.  a=  1.5.  and  nyyo^l.OSo},. 


(6) 

where  do  ^'><1 1’ro  electron  phase  and  axial  velocity  at 
z=0.  The  phase  shift  can  arise  from  the  initial  resonance 
mismatch  (kinematic  shift) 


where  yo  is  the  electron  relativistic  factor  at  z=0  and  the 
changes  in  y  and  u,  result  from  the  electron-wave  interac¬ 
tion  (dynamic  shift) 

S..-("'Ar-AVi.,)£.  (S) 

For  the  moment  consider  only  the  effect  of  the  initial 
mismatch  on  the  electron  azimuthal  phase  bunching.  An 
ensemble  of  electrons  with  initially  equal  transverse  and 
axial  velocities  and  randomly  distributed  in  phase  between 
0  and  2v  if  perturbed  by  an  electromagnetic  wave  which 
satisfies  the  exact  resonance  condition  [Eq.  (S)]  will  bunch 
at  d=ir.  This  is  because  the  electron  initially  at  d=ir  will 
not  be  affected  by  the  wave  while  the  electrons  in  the  neigh¬ 
borhood  of  IT  will  either  slow  down  (slip  in  phase)  if  their 
d  >  IT  or  speed  up  (gain  in  phase)  if  their  d  <  ir  by  the  wave. 
In  this  case  the  electrons  on  the  average  do  not  exchange 
energy  with  the  wave.  Now  suppose  that  Aoj  is  positive. 
From  Eq.  (6)  the  electron  pha%  trajectory  will  have  a 
negative  slope.  This  results  in  forming  the  electron  bunch 
between  0  and  ir  and  the  electrons  on  the  average  convert 
their  kinetic  energy  into  coherent  radiation.  The  efficiency 
of  the  conversion  process  depends  on  where  the  bunch  is 
located.  The  bunch  location  can  be  controlled  by  optimiz¬ 
ing  the  initial  mismatch  (detuning)  of  an  inject  signal  in  an 
amplifier  configuration.’^  It  was  also  demonstrated  that  in 
a  gyro-BWO  configuration’  the  electron  bunch  can  be 
placed  at  d=v/2  by  properly  up-tapering  the  magnetic 
field.  In  doing  so,  the  efficiency  was  increased  from  around 
10%  for  an  untapered  case  to  more  than  30%.  This  was 
accomplished  because  the  up-taper  tends  to  increase  the 
oscillation  frequency  and  the  amount  of  increase  depends 


Table  l.  Pimncicn  used  in  ooinputer 


y 

(beam  cnercy) 

SOkV 

ro 

(tdadvittic  betor) 

1.0M 

/ 

(beam  citimt) 

3  A 

a 

1.5 

(velocity  aptead) 

0. 

(beam  |iiidiii(  center) 

a 

'm 

(waveguide  radius) 

0.113  cm 

on  the  degree  of  tapering.  On  the  other  hand,  in  a  unifmn 
magnetic  field  situation  (gyro-BWO)  the  initial 
may  be  changed  by  varying  the  magnetic  field.  This  arises 
from  the  nonlinear  relationship  between  a  and  when  <o 
is  close  to  the  waveguide  cutoff  frequency.  This  is 
true  as  will  be  verified  in  the  next  section  through  com¬ 
puter  simulations. 

III.  EFFECTS  OF  VARYING  MAGNETIC  RELD  ON 
EFFICIENCY 

One  of  the  advantages  of  utilizing  gyro-BWO  as  a  mil¬ 
limeter  wave  source  is  its  tunability  through  varying  the 
magnetic  field.  Therefore  it  is  important  to  investigate  how 
the  variation  of  magnetic  field  influences  the  output  efS- 
ciency.  In  carrying  out  the  Emulations,  only  the  TEji 
waveguide  mode  was  retained  and  the  parameters  used  are 
shown  in  Table  I  and  the  magnetic  field  was  varied  so  that 
the  electron  cyclotron  frequency  ranged  from  to  20}^ 
where  df  is  the  waveguide  cutoff  frequency. 

An  uniform  section  of  a  waveguide  with  i.=3  cm  was 
employed  and  the  output  power  came  out  from  the  beam 
entrance  plane.  Figure  2(a)  shows  the  output  efficiency 
and  oscillation  frequency  versus  the  imposed  magnetic 
field.  As  is  expected  the  output  frequency  scales  almost 
linearly  with  the  magnetic  field  in  the  entire  region  covered 
in  simulations  except  very  near  the  cutoff.  However,  the 
output  efficiency  exhibits  a  significant  jump  (about  50% ) 
when  the  oscillation  frequency  approaches  the  waveguide 
cutoff  frequency.  Simulation  results  clearly  reveal  that  a 
roughly  constant  output  power  gyro-BWO  can  be  operated 
in  the  following  two  regimes;  high  efficiency  with  about 
8%  tunability  or  medium  efficiency  with  about  40%  tun¬ 
ability.  In  the  remaining  part  of  this  section,  we  shall  elu¬ 
cidate  the  mechanism  which  causes  the  output  efficiency  to 
depend  on  the  electron  cyclotron  frequency. 

Two  examples  will  be  singled  out  to  emphasize  their 
difference.  In  the  simulations  the  electromagnetic  field  is 
decomposed  into  two  independent  components  (Ene  and 
cosine)  in  the  azimuthal  direction.  Figure  3  displays  only 
the  results  of  sine  component  and  the  other  component  has 
similar  behavior.  The  spatial  distribution  of  electric  field  in 
steady  state  for  {l/yQ=  1.6  is  shown  in  Fig.  3(a).  From 
which  the  wave  number  of  the  dominant  mode  is  estimated 
to  be  k(^—0.2li5<i)c.  The  frequency  spectrum  of  the  elec¬ 
tric  field  at  z=:0  is  also  shown  in  Fig.  3(a)  and  is  peaked  at 
<Uo=1.46>c.  Substituting  these  results  into  Eq.  (5)  gives  a 
frequency  mismatch  of  O.OltUf.  By  reducing  magnetic  field 
to  Cle^Yo—  1.05<a,,  the  wave  number  and  frequency  of  the 
dominant  mode  becomes  *^=0-053©^  and  a^s=1.043<Uc. 
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FIG.  2.  Device  petfomunce  venus  mafnecic  field:  (a)  eSkieiicy  i)  (cir> 
cki)  and  oaeiUation  frequency  <uq  (Man)  whicli  ahow  an  incteaae  of  30% 
in  efficiency  when  is  close  to  a,  and  (b)  bunch  phase. 


In  this  case  the  frequency  mismatch  is  0.046d);  which  is 
more  than  four  times  the  previous  case.  This  is  the  conse¬ 
quence  of  nonlinear  dependence  of  o>  vs  (Fig.  1 )  when  a 
is  close  to  Of.  The  effect  of  frequency  mismatch  on  the 
output  efficiency  will  become  clear  after  examining  Fig.  4. 


w 


«/«c 


u/ur 


FIG.  3.  Electric  field  spatial  distribution  and  the  output  frequency  ^wc- 
trum  after  the  system  reaches  steady  state;  (a)  •>■4  (b) 


FIG.  4.  TcM  election  azimuthal  phase  trajectories  after  a>^~2000,  (a) 
n/ro“i-«"r.  (b)  n/ro-»05e),. 


In  order  to  illustrate  how  the  electron  azimuthal 
bunching  process  differs  when  the  magnetic  field  is  varied, 
the  phase  trajectories  of  40  test  electrons  initially  (z— 0) 
uniformly  distributed  between  0  and  2tr  are  followed.  In 
simulations,  the  electron  phases  are  evaluated  at  the  instan¬ 
taneous  position  of  the  electron  and  after  the  system  has 
reached  steady  state.  The  results  of  using  1.6<Uf  are 

given  in  Fig.  4(a).  According  to  £q.  (6)  the  initial  slope  of 
the  phase  trajectory  is  negative  if  the  frequency  mismatch 
( Aa>)  is  positive.  In  the  region  of  rr  <  0  <  2ir,  the  electron- 
wave  interaction  results  in  a  positive  Ay  and  therefore  the 
phase  trajectory  becomes  more  negative  while  in  the  region 
of  0<6<ir  the  negative  slope  should  be  reduced.  These 
detailed  phase  shift  phenomena  are  all  clearly  demon¬ 
strated  in  simulation  results  [Fig.  4(a)].  In  this  case  the 
electrons  eventually  bunch  at  and  zssl  cm.  As  a 
result  of  the  larger  initial  frequency  mismatch  for 
n/yo=1.05(i)(,,  the  electrons  bunch  [Fig.  4(b)]  around 
6=ihr  which  is  substantially  closer  to  d=ir/l  (maximum 
electron-wave  coupling)  than  the  previous  case.  The  vari¬ 
ation  of  the  bunch  phase  versus  magnetic  field  is  shown  in 
Fig.  2(b).  The  resulting  phase  bdiavior  follows  the  effi¬ 
ciency  curve  quite  closely.  This  further  supports  the  asser¬ 
tion  that  the  final  bunch  phase  plays  a  decisive  role  in 
determining  the  output  efficiency  of  a  gyro-BWO.  At  the 
same  time  the  group  velocity  of  the  wave  in  the  latter  case 
is  also  slower  which  increases  the  electron-wave  interac¬ 
tion  time.  These  have  all  contributed  to  the  enhancement 
of  output  efficiency  of  the  latter  case.  An  attempt  has  also 
been  made  to  study  the  effects  of  tapering  the  magnetic 
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field  when  <uo  is  close  to  a),.  Preliminary  results  show  that 
the  output  power  suits  to  oscillate  even  with  a  relatively 
small  up-Upering.  A  more  detailed  investigation  is  needed. 


IV.  SUMMARY 

Computer  simulations  have  been  performed  to  demon¬ 
strate  that  a  gyro-backward  wave  oscillator  operated  at  oig 
close  to  the  waveguide  cutoff  frequency  produces  fu  more 
output  power.  It  was  illustrated  that  one  of  the  plausible 
mechanisms  for  this  high  eflSciency  operation  is  the  result¬ 
ing  substantially  larger  frequency  mismatch  (detuning) 
which  places  the  electron  bunch  close  to  B=ir/2  where  the 
maximum  electron-wave  coupling  takes  place.  To  increase 
the  eflSciency  further,  by  up-Upering  the  magnetic  field, 
more  studies  are  needed. 
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